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1 Why Focus Planet Searches on Binary Stars? 



Searches for planets in close binary systems explore the degree to which stel- 
lar multiplicity inhibits or promotes planet formation. There is a degeneracy 
between planet formation models when only systems with single stars are 
studied — several mechanisms appear to be able to produce such a final re- 
sult. This degeneracy is lifted by searching for planets in binary systems; 
the resulting detections (or evidence of non-existence) of planets in bina- 
ries isolates which models may contribute to how planets form in nature. 
Studying relatively close pairs of stars, where dynamic perturbations are 
the stron gest, provides the mo s t restrictive constraints of this type [see, fo r 
example, iThebault et~aD . booi IPfahl 1200,4 IPfahl and Muterspaughl . I2OO6 . 



In this chapter, we consider observational efforts to detect planetary com- 
panions to binary stars in two types of hierarchical planet-binary configura- 
tions: first "S-type" planets which orbit just one of the stars, with the binary 
period being much longer than the planet's; second, "P-type" or circumbi- 
nary planets, where the planet simultaneously orbits both stars, and the 
plane tary orbital period is much longer than that of the binary Dvorakl . 



1982| . Simulations show each o f these configurations has a large 



of stable configurations [see, e.g. , iBenest 



Pilat-Lohinger and Dvorak*. 2002. Broucke 



20031. iPilat-Lohinger et al 



range 



200 1 , Holr nan and Wiegert 



2003 



1999 



Benestl . ll99& . ,1993> . il989. . ,1988, . iRabl and Dvorak . , 19881 . also, this book, chap- 
ters ZZZZ]. 



2 S-Type Planets 

S-Type planets orbit just one of the stars in a binary, and the binary sep- 
aration is much larger than that between the star and planet. Some of the 
binaries are so widely separated (projected semimajor axis > 1 arcsec- 
ond) that they can be spatially resolved by ground-based telescopes without 
active image correction; for these, traditional planet-finding techniques can 



be used. In fact, astrometric methods often perform best in this regime, as 
the secondary star serves as a convenient reference for the primary, and vice 
versa. Here, astrometric and radial velocity (RV) programs are considered 
as the most versatile search methods. (While transit searches might also be 
possible, these typically have very limited spatial resolutions, and the second 
star can act as a photometric "contaminant.") When the binaries are not 
spatially resolved with simple imaging, modifications must be made to meet 
the measurement precisions required for detecting extrasolar planets. 



2.1 Wide Binaries 

From an observational standpoint, "wide" binaries are considered to be those 
that can be resolved by traditional (uncorrected) imaging techniques. Due 
to atmospheric seeing, this sets the projected sky separation at larger than 
roughly one arcsecond. 



2.1.1 Dualstar Astrometry 



Inter ferometric narrow-angle astrometry phao and Colavital . Il992l . IColavital . 
1994| promises astrometric performance at the 10-100 micro-arcsecond level 
for pairs of stars separated by 1-60 arcseconds. The lower limit of the al- 
lowable binary separation for this technique is that the binary is resolved by 
the individual telescopes in the interferometer; the upper limit is set by the 
scale over which the effects of atmospheric turbulence are correlated. This 
technique w as first de i nonst rated with the Mark III interferometer for short 
integrations Colavital . Il994l |. was extended to longer integrations and shown 
to work at the 100 micro-arcsecond level at the Palomar Testbed Interferom- 
PTI 



eter 



Lane et al.. 2000 



However, achieving such performance requires simultaneous measurement 
of the interferometric fringe positions of both stars, greatly complicating the 
instrument (two beam combiners and metrology throughout the entire array 
are required). In addition, the instrumental baseline vector B connecting 
the unit telescopes must be known to high precision (^ 100 microns). 

In an optical interferometer light is collected at two or more apertures and 
brought to a central location where the beams are combined and a fringe pat- 
tern produced on a detector. For a broadband source of central wavelength 
A and optical bandwidth AA the fringe pattern is limited in extent and ap- 
pears only when the optical paths through the arms of the interferometer 
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are equalized to within a coherence length (A = A^/AA). For a two-aperture 
interferometer, neglecting chromatic dispersion by unequal air paths, the in- 
tensity measured at one of the combined beams is given by 



I{x) 



sin(7rx/A) 
1 + V , i sin (27TX X) 



(1) 



where V is the fringe contrast or "visibility", which can be related to the 
morphology of the source, and x is the optical path difference between arms 
of the interferometer; see Fig. [H More detailed analysis of the operation of 
optical interfer ometers can be found in Principles of Long Baseline Stellar 
Interferometry |Lawsod . |2000| . 



The location of the resulting interference fringes are related to the position 
of the target star and the observing geometry via 

d = B - ^ + Sale's, +c (2) 

where d is the optical path-length one must introduce between the two arms 
of the interferometer to find fringes (often called the "delay"), S is the 
unit vector in the source direction, and c is a constant additional scalar 

delay introduced by the instrument. The term 6a(^S is related to the 

differential amount of path introduced by the atmosphere over each telescope 
due to variations in refractive index. 

If the other quantities are known or small, measurement of the instru- 
mental path length d required to observe fringes determines the position of 
the star S . For a 100-m baseline interferometer, an astrometric precision 
of 10 /ias corresponds to knowing d to 5 nm, a difficult but not impossible 
proposition for all terms except that related to the atmospheric delay. Atmo- 
spheric turbulence, which changes over distances of tens of centimeters and 
on millisecond timescales, forces one to use very short exposures to maintain 
fringe contrast, and hence limits the sensitivity of the instrument. It also 
severely limits the astrometric accuracy of a simple interferometer, at least 
over large sky-angles. 

However, in narrow-angle astrometry one is concerned with a close pair of 
stars, and the observable is a differential astrometric measurement, i.e. one 
is interested in knowing the angle between the two stars (A^ = 52 — si)- 
The atmospheric turbulence is correlated over small angles. If the measure- 
ments of the two stars are simultaneous, or nearly so, the atmospheric term 
subtracts out making possible high precision "narrow-angle" astrometry. 
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Figure 1: The response of an interferometer. The top two curves have been 
offset by 2 and 4 for clarity. The widths of the fringe packets are determined 
by the bandpass of the instrument, and the wavelength of fringes by an 

averaged wavelength of starlight. The top curve shows the intensity pattern 
obtained by observing two stars separated by a small angle on the sky — the 
observable is the distance between the fringe packets. 



4 




common metrology fiducials 



h > 


Beam Combiner 1 


L A 






V * 


* vl 



Held separator 



star selector 
mirror 



Beam Combiner 2 




Figure 2: Schematic of splitting the hght in a dualstar interferometer. 



The requirement that the target and reference stars be observed simul- 
taneously results in a significant instrumental complexity, i.e. essentially 
two complete interferometers are required to share the same set of apertures 
(see Fig. [2]). The splitting of light from the stars into two separate sets of 
delay lines, beam transport systems and beam combiners is done in a "dual- 
star module" located just after the apertures, with the split generally being 
accomplished using a beam-splitter. Considerable care must be taken in 
designing the system in order to avoid small pathlength measurement errors. 

The exact level of astrometric precision that can be achieved depends on 
many factors, including the separation of the target /reference pair, the size 
of the interferometric baseline and the levels and distribution of atmospheric 
turbulence. For a typical Mauna Kea seeing profile the astrometric precision 
is 

9 

— 300^ — - arcsec (3) 

where B is the baseline length in meters, 9 is the target /reference separation 
in radians, and t is the integration time in seconds. For typical baselines of 
~ 100 m, and an angular separation of ~ 30 arcsecond implies an astrometric 
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Figure 3: Astrometric accuracy vs. star separation in a one-hour integration 
for different baseline lengths. Model atmospheres providing 1/2- and 1.0- 
arcsecond seeing are shown. These results assume an infinite outer scale, 
and better results are achieved when the baseline exceeds the outer scale, as 
would be expected with a 100 m baseline at most sites. Measurements with 
the Mark III interferometer of a 3.3 arcsecond binary sta r are consistent with 
the model. This figure is from IShao and Colavital 1992 . 



precision of 30 /las in an hour (see Fig. [3]). 

The magnitude of the astrometric signal of the star's motion about the 
center of mass (CM) between it and its planet is given by: 

AacM = 2^a, = »^. (4) 
^ Ms M,/M0 524 ^ ' 

where Mp, M;,, Mj, M© are, respectively, masses of the planet, star, Jupiter, 
and Sun, and ap is the semimajor axis of the planet's orbit. 
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Thus, the minimum mass that can be detected is roughly 



Mp/Mj 



> 524- ' 



> 



Mp/Me > 1.6 



^ ^ cra/20/ias d Ms 
ap/lAU lOpcM© 
cTa/l/xas d Ms 



ap/lAf/lOpcM 







(5) 
(6) 
(7) 



where M^ is the mass of the Earth, and here d is the distance to the target 
star. 



2.1.2 Radial Velocities 

When the stars in a binary can be spatially resolved without active image cor- 
rection on ground based telescopes, the spectrum of each star can be recorded 
separately without contamination from the other, and the st andard precision 



RV method described be low can be used [see, for example, ICampbell et al. 



19881 . iButler et al.l . Il996| . Similarly, if the secondary is much fainter than 
the primary, precision RV might be performed on the brighter star as though 
it were single, though there is concern about the influence of the fainter 
lines. Several (~ 30)exoplanet candidates in binaries have been discovered 
in this manner. In some cases the stars were not previously known to be 
binaries, and their natures were only discovered by long-term RV trends or 
follow-up adaptive optics imaging. So me of these efforts to detect planets 
in binary stellar systems i nclude that of Tovota et al.l 2005| for single-lined 



and wide binaries, that of Desidera et al. 



20061 targetir ig wide binaries, a nd 



the program targeting single- lined and wide binaries of lUdry et al.l [2004 . 

The highest precision RV observations are obtained either from the I2 
(molecular iodine) absorption cell or the use of carefully designed spectro- 
graphs with fiber scrambling. In order to achieve a RV precision of ~ Im 
an iodine absorption cell is used to superimpose a reference spectrum on the 
stellar spectrum (by sending a starlight through the cell). The spectrum 
provides a fiducial wavelength scale against which radial velocity shifts are 
measured. 

Thanks to its conceptual simplicity, the iodine technique is the most com- 
monly adopted way to obtain precision radial velocities. Iodine absorption 
cells are available on many spectrographs — HIRES at the 10m Keck I (Keck 
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Observatory), Hamilton at the 3m Shane (Lick Observatory), SARG at the 
3.6m TNG (Canary Islands), UCLES at the 3.9m Anglo- Australian Telescope 
(Anglo-Australian Observatory), HRS at the 9m HET (McDonald Observa- 
tory), MIKE at the 6.5m Magellan (Las Campanas Observatory), UVES at 
the 8m Kueyen (Cerro Paranal), HDS at the 8.2m Subaru (National Astro- 
nomical Observatory of Japan) and many other — and are used for planet 
detections. 

In the iodine absorption cell technique, the Doppler shift of a star spec 



trum is determined by solving the following equation [Marcy and Butlerl . 



1992 



lobsW = [IsiX + AA,) Tj.iX + AA,J] ® PSF (8) 

where AA^ is the shift of the star spectrum, AXj^ is the shift of the iodine 
transmission function Tj,^, represents a convolution, and PSF a spectro- 
graph's point-spread function. The parameters AA^, AA72 well as param- 
eters describing the PSF are determined by performing a least-squares fit to 
the observed spectrum, lobs, as seen through the iodine cell. To this end, 
one also needs a high SNR star spectrum taken without the cell, Ig, which 
serves as a template for all the spectra observed through the cell, as well 
as the I2 transmission function, Tj^, obtained with the Fourier Transform 
Spectrometer at the Kitt Peak National Observatory. The Doppler shift of 
a star spectrum is then given by AA = AAs — AXj^. 

The velocity reflex amplitude of a star due to an unseen companion is 
given by 

2'Ka„smiri M„ 2v^M„sini„ 
Avb = 2-^^ P - VP 



Pp Ms + Mp ^{Ms + Mj 



pj U/p 



-1 (Mpsimp/Mj) 
= 56.9ms X , (9) 
^{(Ms + Mp) /Mo)(ap/lAU) 

where Pp is the period of the planet's orbit, G is the gravitational constant, 
and ip is the inclination of the planet's orbit. 

For an RV precision of arv, the minimum mass that can be detected is 
roughly 

Mpsinip/Mj > 0.018^-^y/ ((M, + Mp) /Mq) (a^/lAU) (10) 
Mpsin^p/Me > 5.6^^^^((M, + Mp) /Mq) (a^/lAU). (11) 
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Figure 4: Sensitivity to S-Type planets in wide binaries, comparing astromet- 
ric and radial velocity techniques. All calculations assume solar mass stars; 
astrometric sensitivity assumes a distance of 10 pc to the target system. 

2.1.3 Observational Precisions 

Astrometry is most sensitive to long period planets, RV to short period ones. 
Figure H] shows the companion masses one can detect for each method, as- 
suming 20 /ias ground-based astrometry, 1 fias space-based astrometry, and 
lms~^ RV precisions. 

2.2 Close Binaries 

Radial velocity surveys for extrasolar planets have been restricted largely to 
stars within ~100 pc, and the majority of detected exoplanets are at distances 
of 10-50 pc. Therefore, our observational definition of a "wide" binary cor- 
responds to projected orbital separations of >50AU. More compact systems 
form the complementary class of "close" binaries. From a theoretical stand- 
point, binaries with semimajor axes of <50AU pose important challenges to 
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standard ideas about the formation of giant planets. 

Imagine a protoplanetary disk around one star in a newly formed binary. 
Suppose that the orbit has semimajor axis a and eccentricity e, and assume, 
for simplicity, that the stars have equal masses. The tidal gravitational field 
of th e companion star trun cates the disk at a radius of Rt — 0.26a(l — e^)^'^ 



Pichardo et al.. 2005 . If Rf < 3 AU (i.e., inside the so-called "ice line'' 



[e.g. 

it seems unlikely that icy grains could form and grow into planetesimals, thus 
precluding the embyoni c stage of giant planet formation in the core-accretion 



scenario [e.g., iLissauerl . Il993| . In more extreme cases, when Rt < 1 AU, 



there may be i nsuffi cient material in the disk to yield a Jovian-mass planet 



Jang-Condelll . l2006l |. Even when Rt is as large as ~10AU, stirring of the 
disk by the tidal field and the thermal dissipation of spiral waves may inhibit 
lanetes i mal formation, as wel l as stabilize the disk against fragmentation 

In this case 



i )ianetes i mai lormation, as wel l as si 
Nelsonl . l2000l . Ixhebault etaP . 12004 



2006 



it may b e that 

neith er the core-accretion picture nor gravitational instability [e.g., iBosd . 
2000| are accessible modes of giant planet formation. We adopt Rt = 10 AU 
(a < 40 AU for modest e) as a fiducial upper limit for which Jovian planet 
formation is significantly perturbed and perhaps strongly inhibited. This 
provides a simple theoretical definition of a close binary that roughly matches 
our observational measure. 

A handful of planets in binaries with a < 20 AU have already been dis- 
covere d (see Tabled]). The tightest of these systems [HD 188753 iKonackil . 
2005bl has a periastron separation of only ~6AU, which seems severely at 



odds with the conventional lore on Jo vian planet formation. iPfahll 2005| and 



Portegies Zwart and McMillanI |2005[ | suggested that the planetary host star 



may have been acquired in a dynamical exchnage interaction in a star cluster 
after the planet formed, thus circumventing the complicating factors listed 
above. As most stars are born in clustered environments, one wonders how 
often dyn amics can account for p l anets in close binaries. This idea was ex- 
plored in iPfahl and Muterspaughl |2006l | , where it was found that exchange 
interactions can account for only ~0.1% of close binaries hosting planets. 
However, the (admittedly small) sample of systems in Table 1 seems to in- 
dicate that a larger fraction of ~1% of close binaries harbor giant planets 
(see Pfahl & Muterspaugh 2006 for details), and perhaps planets do some- 
how form frequently in these hostile environments. It crucial that we begin 
to develop a census of planets in close binaries in order to test the different 
theories about planet formation and dynamics. 
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Close Binaries with Planets 


Object 


ab (AU) 


e Ml /Ma Rt 


(AU) Refs 


HD 188753 


12.3 


0.50 1.06/1.63 


1.3 1 


7 Cephei 


18.5 


0.36 1.59/0.34 


3.6 2, 3 


GJ 86 


~20 


0.7/1.0 


~5 4, 5, 6 


HD 41004 


~20 


0.7/0.4 


~6 7 


HD 196885 


~25 


1.3/0.6 


~7 8 



Table 1: When no eccentricity is given, only the projected binary separation 
is known. M1/M2 is the planetary host mass divided by companion mass. In 
HD 188753, the secondary is a binary with semimajor axis 0.67 AU. In GJ 
86, the secondary is a white dwarf; to estimate the ti dal truncation ra dius 
Rf , an original compan ion mass of IMf^ is assumed. flllKonackil f2005bl: (2) 



Campbell et al.l 11988^ f3l Hatzes et all 120031: f4l Queloz et all 120001 (5) 



Mugrauer and Neuhauser 2005] (6) Lagrange et al. 2006| (7)~ Zucker et al 



2004 f8) IChauvin et al.l |2006l 



2.2.1 PHASES Astrometry 



The dualstar astrometry method can be modified for use when the binaries 
are so clo se that the individual tel e scope s of an interferometer cannot resolve 
the pair Lane and Muterspaughl . |2004| . The interferometer itself overre- 
solves the binary, as in Fig. [T| its high spatial resolution then allows for 
precision astrometric measurements. 

In this mode, the small separation of the binary results in both compo- 
nents being in the field of view of a single interferometric beam combiner. 
The fringe positions are measured by modulating the instrumental delay with 
an amplitude large enough to record both fringe packets. 

However, since the fringe position measurement of the two stars is no 
longer truly simultaneous it is possible for the atmosphere to introduce path- 
length changes (and hence positional error) in the time between measure- 
ments of the separate fringes. To reduce this effect, a fraction of the incom- 
ing starlight is redirected to a separa te beam- combiner . This beam-combiner 
is use d in a "fringe-tracking" mode Shao and Staelin . 1980l . Colavita et al. 



1999 



where it rapidly (10 ms) measures the phase of one of the starlight 
fringes, and adjusts the internal delay to keep that phase constant. The 
fringe tracking data is used both in real-time as a feed-back servo, after 
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Figure 5: The expected narrow- angle astrometric performance in mas for 
the phase-referenced fringe-scanning approach, for a fixed delay sweep rate, 
and an interferometric baseline of 110 m. Also shown is the magnitude of 
the temporal loss of coherence effect in the absence of phase referencing, 
illustrating why stabilizing the fringe via phase referencing is necessary. 

which a small residual phase error remains, and in post-processing where 
the measured residual error is applied to the data as a feed-forward servo. 
This technique — known as phase referencing — has the effect of stabilizing the 
fringe measured by the astrometric beam- combiner. For this observing mode, 
laser metrology is only required between the two beam combiners through 
the location of the light split (which occurs after the optical delay has been 
introduced), rather than throughout the entire array. Without phase refer- 
encing, the astrometric precision obtainable is a factor of a hundred worse; 
see Fig. [51 

To analyze the data, a double fringe packet based on Eq. [T] is then fit to 
the data, and the differential optical path between fringe packets is measured. 
A grid in differential right ascension and declination over which to search is 
constructed. For each point in the search grid the expected differential delay 
is calculated based on the interferometer location, baseline geometry, and 
time of observation for each scan. A model of a double-fringe packet is then 
calculated and compared to the observed scan to derive a value; this is 
repeated for each scan, co-adding all of the values associated with that 
point in the search grid. The final surface as a function of differential R.A. 
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and declination is thus derived. The best-fit astrometric position is found at 
the minimum-x^ position, with uncertainties defined by the appropriate 
contour — which depends on the number of degrees of freedom in the problem 
and the value of the x^-minimum. The final product is a measurement of the 
apparent vector between the stars and associated uncertainty ellipse. Because 
the data were obtained with a single-baseline instrument, the resulting error 
contours are very elliptical, with aspect ratios that sometimes exceed 10:1. 

The Palomar High-precision Astrometric Search for Exoplanet Systems 
(PHASES) program uses this technique to monitor ~ 50 binaries to search 
for substellar companions, k Pegasi is a well-known, nearby triple star sys- 
tem. It consists of a "wide" pair with semi-major axis 235 mas (8.14 AU), 
one component of which is a single-line spectroscopic binary (semi-major 
axis 2.5 mas, physical separation 0.087 AU\ Fig. [6]). The perturbation due 
to the unseen (faint) short-period component is evident; similar sized pertur- 
bations with longer orbital periods would indicate the presence of planetary 
companions. Figure [7] shows the mass-period phase space in which PHASES 
observations show companions do not exist in face-on, circular orbits in the 
13 Pegasi system. 



Figure 6: The visual orbits of k Pegasi showing perturbations by the third 
component. The spiral line represents the apparent motion of the short- 
period pair's center of light; the ellipses represent the la uncertainties for 
PHASES measurements. 
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Figure 7: The 13 Pegasi Mass- Period companion phase space shows PHASES 
observations can rule out tertiary objects as small as two Jupiter masses. A 
few mass-period combinations introduce slight improvements over the single- 
Keplerian model, but none of these are more significant than 1.7a, and are 
probably not astrophysical in origin. There is a long period cutoff in sen- 
sitivity due to the finite span of the observations. Similar detection limits 



have recently been published on other PHASES targets [Muterspaugh et al. 
20061. 



2.2.2 Radial Velocities 

In the case when a composite spectrum of a binary star is observed, the 
classical approach with the iodine cell (described in section I2.1.2p cannot be 
used since it is not possible to observationally obtain two separate template 
spectra of the binary components. To resolve this problem, one can proceed 
as follows. First, one always takes two sequential exposures of each (binary) 
target — one with and the other without the cell. This is contrary to the 
standard approach for single stars where an exposure without the cell (a 
template) is taken only once. This way one obtains an instantaneous template 
that is used to model only the adjacent exposure taken with the cell. Next, 
one performs the usual least-squares fit and obtain the parameters described 
in eq. [Hi Obviously, the derived Doppler shift, AAj (where i denotes the epoch 
of the observation), carries no meaning since each time a different template is 
used. Moreover, it describes a Doppler "shift" of a composed spectrum that 
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is typically different at each epoch. However, the parameters — in particular 
the wavelength solution and the parameters describing PSF — are accurately 
determined and can be used to extract the star spectrum, /*^*(A), for each 
epoch i, by inverting the eq. [8j 



'12) 



where denotes deconvolution, and PSF^ represents the set of parame- 
ters describing PSF at the epoch i. Such a star spectrum has an accurate 
wavelength solution, is free of the I2 lines and the influence of a varying PSF. 
In the final step, the velocities of both components of a binary target can 
be measure d with the well know two- dimensional cross- correlation technique 
TODCOR ^Zucker and Mazehl . Il994l | using as templates the synthetic s pec- 
tra derived with the ATLAS 9 and ATLAS 12 programs Kuruca . Il995| and 
matched to the observed spectrum, /^(A). The formal errors of the velocities 
can be derived from the scatter between the velocities from diff e rent e chelle 
orders or using the formalism of TODCOR Zucker and Mazehl . Il994| . The 
technique c urrently produces RVs of binary stars with an average precision 



of 20 ms [Konackil . l2005a| . Improvements to the technique are being in- 
troduced to reach the level below 10 ms~^. 



With the iodine technique iKonackil |2005al | has initiated the first RV sur- 
vey for circumprimary or circumsecondary planets of binary or multiple stars 
(mainly hierarchical triples). The survey's sample of ~ 450 binaries (north- 
ern and southern hemisphere) was selected based on the following criteria. 
(1) The apparent separation of the components had to be smaller than the 
width of the slit (e.g., 0.6 arcseconds for Keck-I/HIRES) to avoid possible 
systematic effects. Such systems will remain unresolved under most seeing 
conditions. (2) The brightness ratio between the components should not be 
too large (at the order of 10 or less) to be able to clearly identify spectra 
of both components. (3) The orbits of the binaries should be well known 
to constitute a firm ground on which one can discuss possible detection (or 
lack) of planets and substellar companions in the context of the binary char- 
acteristics. All these requirements could be satisfied by targeting a subset of 
speckle binaries th at have determined orb its from the Catalog of Orbits of Vi- 



sual Binary Stars [Hartkopf et al.l . l200l| . containing 1700 binaries, of which 
1300 have projected semi-major axes smaller than 1 arcsecond. The sample 
is sufficiently large to produce meaningful statistics. Also, such binaries have 
been ignored by previous RV studies. 
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The survey was initiated in 2003 at the Keck-I/HIRES and so far has 
covered about 50% of the northern hemisphere target sample, of which ~150 
stars have been observed at least twice. In July 2005, t he first candidat e 



planet in a triple star system HD 188753 was announced [Konackil . 12005b 
A few more candidates as well as HD 188753 are currently being scrutinized. 
In the course of the Keck survey, a dozen of new triple star systems with the 
third bodies in the low-mass star or even brown dwarf regime have been also 
identified and will be published in the near future. 

2.2.3 Eclipse Timing 

Should a binary happen to be oriented with its orbital plane in the line of 
sight, it will exhibit eclipses as one star passes in front of another. For a 
binary with separation wide enough to allow for stable planetary systems 
to exist around just one component, the probability of such an alignment is 
extremely small, and very few targets are accessible. However, should one 
find such a fortunate happenstance, one can detect the planetary companions 
by precision timing of the eclipses. Clearly, this method has no direct analog 
for single systems. 

This method can detect S-type planets or similarly moons of transiting 
planets. For this evaluation it is assumed that the depth of the planet (or 
moon) eclipse is sufficiently small as to be ignored (in such a detection, one 
could then reevaluate light curves to look for such transit signals) and the 
binary orbit is circular. The velocity of the binary orbit and the offset of the 
star-planet CM from that of the star by itself determines the timing variation 
observed as 

At = XcM/Vb 

' ttp sin (pMp \ f Pb 



M2 + Mp J \2nab 

, , , N Op iM^/Mj) , , 

^ 57 seconds x P^/month /»w,/,,. ^ sin0 13 

(ab/7) (M2/M0) 

where M2 is the mass of star 2 (or the transiting planet, assumed to host 
the S-type companion), Mp is the mass of the S-type object orbiting M2, 
is angle between the planet's orbital angular momentum vector and the 
direction of motion of the host star during eclipse, and Pb is the period of 
the binary orbit. The timing delays due to orbit of M2 about the M2-Mp 
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CM. The factor of 7 appears in the final form is an approximate criteria 
for stabihty, that the planet semimajor axis is 7 times smaller than that of 
the binary (this factor varies by the system, and can be determined through 
detailed simulations); the above is thus an upper limit for the timing effect. 
Converting the se mimajor axis to orbital periods, 

At 41 seconds x (Pb/month) ^ (P^/day) ^ - — ^^^^'^ sin 



1 2 Mp/Mj 

65 seconds x {Ph / month) (Pp/day)^ M /M. 



(14) 



where — Mi + M2 + Mp. The final form assumes Mi M2, in which case 
the maximum stable planet period is a thirteenth that of the binary period, 
implying days and months are the natural units for each respectively (S- 
type planets cannot exist in much shorter period systems, and longer period 
systems are even less likely to show eclipses). The equivalent relationship for 
a moon orbiting an eclipsing Jupiter is 

At ^ 13.3seconds x (Pfc/month)^ (P /day)^ Mp/M^^ -sincf) 

(M6/Mo)3 (M2/Mj)3 

(15) 

where now the h subscript refers to the star-Jupiter analog system and p to 
the Jupiter analog's moon. 

The precision with which eclipse minima can be timed is derived using 
standard fitting techniques. Assume a photometric data set {yj} occurring 
at times {ti} with measurement precisions {cTj}, and a model photometric 
fight curve of fiux F{t — to). The corresponding intensity is I{t — to) — 
fF{t — to)7TD'^At/4, where / (0 < / < 1) is the fractional efficiency and 
throughput of the telescope, D is the telescope diameter, and At is the 
sample integration time. (P(t — to) might be determined to high precision by 
observing multiple eclipse events.) The fit parameter to is uncertain by an 
amount equal to the difference between the value for which is minimized 
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and that for which it is increased by one: 1 + x^ito) = X^i^o + (^to), 



N 



1=1 



Vi - I{ti- to) 



N r 



N 

E 

i=l 



Vi - I{ti- to) 



1 2 



+ 



E 

i=l 

N 

E 

8=1 



dl(t) , 



Vi - I {ti - to - ato) 



yi-I {ti - to) - ( 



dlit) 



-1 2 



(^). . . [y^-I{t^-to)]at, 



Because to is the minimizing point, the first derivative of at to is zero, 
giving 



dx^ it) 
dt 



t=to i=l 

Rearrangement of terms leads to 

1 



x^(^). . . iy^-Ht^-to))at,, 
2^- ^*=*'"*° = 0. (16) 



'E 



N 



N ^ fdmy 

/ t=ti—to 



An echpse of length r is approximated as a trapezoid-shape light curve 
(see Figure[8]) with maximum and minimum photon fluxes Fo and Fo{l — h/2) 
[h is a dimensionless positive number producing an eclipse depth of hFo/2; 
in the case of a faint secondary, h is roughly twice the ratio of the squares of 
the stellar radii, 2R2/RI). The ingress and egress are each assumed to be of 
length kT/2 [k ^ 2i?2/(-Ri + R2) is unity in the case of an eclipsing binary 
with components of equal size, when the trapezoid becomes a "V" -shape). 
In functional form, this model is: 



r Fo 



F(t-to) 



Fo{l-ht/{kT)-h/{2k)) 
Fo{l-h/2) 

Fo{l + ht/{kT)-h/ {2k)) 



t - to < -r/2 
-r/2 < t-to < -T/2 + kT/2 
-T/2 + kT/2< t-to <T/2-kT/2 
T/2-kT/2< t-to <r/2 
r/2 < t - to 

(17) 
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Model lightcurve for an eclipsing binary 




-x/2 -(1-k)T/2 (1-k)T/2 t/2 

Figure 8: Eclipsing binary model light curve. 

(also, see Figure [H]). 

Only the portions of the light curve during ingress and egress have nonzero 
(in more accurate models, the light curve slope will be nonzero but 
small in other regions, and will not contribute much to the sum); this slope 

is\^\ = hFo/{kT). 

The number of data points contributing to the sum is thus N = gkr/At, 
where < 5^ < 1 is the fraction of the eclipse observed (and also accounts for 
the fraction of time lost to, e.g., camera readout) and At is the integration 
time for each measurement. 

The measurement noise ai is given by 

crj= {1 + al^ + hg + HdarkAt + al^) ^ (18) 

where hg = fFbgiiD^At/A is is the sky background, Udark is detecter dark 
current, arn, is detector read noise, and asc is scintillation noise given by 



YoungI |1967 



as 



asc = 0.09/ (D/lcm)"^/^Xe~'^/(^°°°'^V(2At/l second)^ (19) 

^ 0.003/ (D/lm)"^/V(AVl second)^ (20) 

where X is the airmass and h is the altitude of the observatory. The drop in 
noise during eclipse is ignored (a factor less than ^ 1.4) and equation [TS] is 
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combined with equation [T7] to obtain an overall timing precision (in seconds) 
of 



fc (r/l second) / 4 9x10"'' wD^ fF^g/i + n^ark + second)" 



^ 0.18 seconds X J o o" H — J7¥ ■ (21) 

The final term in the first line (associated with dark current, read noise, and 
background) is generally smallest and will be ignored. In most cases, the 
second term — associated with scintillation — is dominant (though zero in the 
case of space-based observatories). The exponent of {V — 12)/2.5 shows that 
for meter-sized telescopes, photon noise is only dominant for stars fainter 
than twelfth magnitude. 

Systematic and astrophysical noise sources may have effects that limit 
the actual precisions achieved. Mass transfer between stars can cause drifts 
in orbital periods. Variations of this type are non-periodic, distinguishing 



them from companion signals. lApplegatd 1992l | has shown that gravitational 



coupling to the shapes of magnetically active stars can cause periodic modu- 
lations over decade timescales. This mechanism requires the star to be inher- 
ently variable; false positives can be removed using the overall calibrations of 
photometric data. It is possible that star spots will hav e large effects on tim- 



ing r esiduals that are particularly difficult to calibrate [Watson and Dhillonl . 



2004| . Due to orbit-rotation tidal locking, the effect of a starspot on the light 
curve can be detected from the light curves of several orbits, and starspot 
fitting potentially can remove the timing biases introduced. 

The timescale for eclipses of such long period binaries is of the order of 12 
hours. Comparing eqs. [H] and [211 one finds that S-type planets with periods 
of a few days can be detected around either star in eclipsing binaries with 
month-long periods if they mass more than 



(Pfo/month)^ (Pp/day)^ V \ (D/lmf {D/lmf 

(22) ' 

Eq. [21] indicates that a meter-class ground-based telescope can time a 
giant planet transit {h ^ 0.02, k ^ 0.18) to approximately 9.4 seconds in the 
regime where the photometric precision is dominated by scintillation noise, 
assuming a Jupiter sized planet orbiting a star of solar size and mass with 
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period of a month (implying 6-hour duration echpses). This precision is suf- 
ficient to find Earth mass moons. For bright stars, space-based observatories 
offer even better precisions. Unfortunately, no transiting exoplanets with 
periods this long have yet been discovered. 

Space-based photometric missions such as Kepler have as their primary 
goal the detection of Earth-like planets via transits of the planet across the 
star. However, such photometric events can be explained by other astrophys- 
ical phenomena, such as a transiting Jupiter blended with a background star, 
so these results may be unreliable. However, Earth-like moons of transiting 
Jupiters might be identified through timing, and it is possible to confirm 
the nature of such a system. In such a scenario, a transiting Jupiter can 
be positively confirmed by ground-based radial velocity observations. Once 
this has been established, variations in the transit times would be used to 
detect Earth-sized moons. Because these photometric missions have limited 
lifetimes (~ 3 years), detections of moons are only possible for short period 
(few months or less) Jupiters, for which many transit events can be observed 
(unless a follow-up ground-based campaign is pursued with large telescopes). 
If the planet/moon are to be in the habitable zone, one must look for such 
systems around late-type (cool) stars. It is possible that such systems have 
the greatest likelihood of being habitable; tidal-locking of the Earth-sized 
moon to the Jupiter-like planet would ensure that the moon has day/night 
cycles and stabilize its rotational axis similar to the way in which the Earth's 
is stabilized by its own moon. Bot h of these condition s have been argued as 



favorable for life [see, for example, iLaskar et al.l . 11993 



2.2.4 Observational Precisions 

Figure [9] shows the companion masses one can detect for each method, as- 
suming 20 fias ground-based astrometry, 20m s~^ RV, and 1 m ground-based 
photometric telescope for eclipse timing. 



3 P-Type (Circumbinary) Planets 

All the confirmed planets found in binary systems thus far are in S-type 
orbits. Discovery of circumbinary planets would consitute a new class of 
solar system, and would inspire new considerations to the interplay between 
system dynamics and planet formation. 
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Figure 9: Sensitivity to S-Type planets in narrow binaries, comparing astro- 
metric, radial velocity, and eclipse timing techniques. All calculations assume 
solar mass stars. The PHASES sensitivity assumes 20 /xas precision and a 
distance to the system of 20 pc. Eclipse timing assumes aim photometric 
telescope observing a, V = 10 magnitude system with binary orbital period 
2 months (longer period systems are even less likely to show eclipses); the 
eclipse timing sensitivity curve only extends to the region where planets are 
likely to have stable orbits. 
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3.1 Radial Velocities 



A circumbinary planet will exhibit two indirect effects on the velocities of the 
stellar components of the system. First, the apparent system velocity will 
vary in a periodic manner due to the motion of the binary about the system 
barycenter. Second, the finite speed of light will cause apparent changes in 
the phase of the binary orbit. These effects may be detectable using modern 
observational techniques. 

The first effect is that the binary will exhibit periodic changes in the 
apparent system velocity; this is the same effect as seen in a single star. 
However, it may be harder to detect for three reasons: (1) the binary system 
is usually more massive than a single star of the same magnitude, (2) ex- 
tremely short-period planet orbits (to which system velocity measurements 
are most sensitive) are unstable around binaries, and (3) the presence of two 
sets of spectral lines may complicate the measurement, as in section 12.2.21 
Equation [9] shows that a Jupiter massed planet with the shortest period sta- 
ble orbit around a 10-day period binary causes a 2Mq binary to move about 
its barycenter by ~ 40ms~^, with the amplitude decreasing as the square 
root of planet orbit semimajor axis. Radial velocity observations with the 
20ms~^ precision demonstrated with Konacki's method can detect Jupiter- 
like planets in orbits of size ~ 4 AU or less, down to the critically stable 
orbit. 

The second observable effect is the additional light travel time as the 
binary system undergoes reflex motion caused by the planet. The magnitude 
of this effect is given by 

At = 2^^^^^^^ = 0.95 seconds x K/1 AU) (M,/M,) sin 
cMfe Mb/M0 

Following a similar derivation as that for finding the expected precision 
of eclipse timing, one finds the precision with which one can estimate the 
orbital phase of a binary based on radial velocity measurements is 

(23) 



where arv is the radial velocity measurement precision and ^ is the deriva- 
tive of the model radial velocity curve with respect to orbital phase, evaluated 
at times tj. The timing precision corresponding to the phase precision derived 
is given by ^ = ^. 
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Approximating the binary orbit as circular, v{t) K cos + (fj- K 

N measurements (each with two measured velocities, one for each star) are 
approximately evenly distributed in phase, 

^ (24) 

^ PWrv ^25) 

^2 (2Ar-12)7rX ^ ' 

where 12 is the number of degrees of freedom for the model. 

If the lines from both stars are observed, the effective K ys, K2 and 
the resulting (Icr) minimum detectable mass is thus 



(a^^/20ms-i) (PjlOdays)^/^(Mb/M0)^/^ , , 

M„ = 41.4Mj X ^ ' , — ^ ' ^' — , 26 



where ij, and ip are the inclinations of the binary and planet orbits, respec- 
tively. Twenty- five 20ms~^ radial velocity measurements of the "prototypi- 
cal" system could detect moderate-mass brown dwarfs (ps 30Mj) at critical 
orbit. Objects at the planet/brown dwarf threshold of 13 Mj are only de- 
tectable in orbits larger than 0.82 AU around a ten-day binary of sunlike 
stars. Alternatively, if only one set of lines are observed, the resulting ex- 
pression is 

/ Mi\ (c7rJ20ms-i) (P6/lOdays)^/^(M6/M0)2/^ , , 

M„ = 41.4Mj X 1 + — P^^^^^^ — ' ^ ' — ^ ' ^' — , (27) 

^ \ M2J ViV - 11 sin % sin ip (op/1 AU) ^ ' 

where Mi is the mass of the star whose lines are observed, and M2 is that 
of the faint star. 

High precision radial velocity observations are only possible on slowly 
rotating {vsini < 10ms~^) stars; measurements of more rapidly rotating 
stars are limited by line broadening to levels worse than the nominal 20 m s~^ 
that has been referenced by this work. This effect is particularly important 
for finding planets around short-period binaries, in which the stars' rotation 
rates are often tidally locked to the binary orbital period; these rotation 
rates limit the observed precisions for systems with periods approximately 
five days or less. 
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Figure 10: Sensitivity of radial velocity measurements to circumbinary plan- 
ets. The two vertical lines at the left represent the approximate critical orbits 

around 5-day (to the left) and 10-day period binaries; shorter period com- 
panions have unstable orbits. Stars whose rotation rates are tidally locked to 
orbital periods less than about 5 days show sufficient rotational line broad- 
ening to prevent 20ms~^ radial velocity precisions. The calculations assume 
the binary consists of two stars each massing 1 M©. 
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3.2 Eclipse Timing 



It has long been recognized that periodic shifts in the observed times of 
photometric minima of eclipsing binaries can indicate the presence of an ad- 
dition a l component to the system ( see, for example. IWolti e3 |l922t . Ilrwini 
1952], Frieboes-Conde and Hercz"e3 |l973 |. Doyle et al. 199^). The ampli- 
tude of the effect is given by eq. [231 As with RV measurements, there is a 



mass/inclination ambiguity; the following derivation assumes no correlation 
between binary and planet inclinations. 

Dividing the precision of an individual measurement by Nobs — 6 (where 
Nobs is the number of eclipses observed and there are 6 parameters to a 
timing perturbation fit, two periods and the eccentricity, angle of periastron, 
epoch of periastron, and mass ratio of the wide pair), converting Fq to V 
magnitude, and combining eq. [21] with that for the timing effect of reflex 
motion (eq. [231) gives a minimum detectable companion mass of 



Mp = 0.19MjX. 



k (r/1 hr) 



Mb/Me 



;Lo(V-12)/2.5 



/ 



f gh^Nobs - Q) {ap/1 AV) sin ip \ (D/lm) 



(28) 

One might also inquire about the sensitivity of this technique to outer 
planets in systems comprised of a single star and a transiting "hot" Jupiter. 
In this case, h ^ 2Rl/R%^^ ^ 0.02 and k ^ 2Rp/{R,tar + Rp) ~ 0.18, the 
"binary" is half as massive, and the eclipse duration is half as long. The 
companion sensitivity drops by a factor of 8, and the technique is (barely) in 
the range of detecting additional companions of planet mass. However, for 
the typically V = magnitude transiting planet systems being discovered, 
3 ms~^ radial velocity observations are more sensitive than half-meter tele- 
scope transit timing for companions with periods up to 60 years; even for 
observatories such as HST and SOFIA (for which scintillation noise is small 
or zero), this transition occurs at 15 year period companions. 

It should be noted that the above description does not account for the 
possibility of resonant orbits, for which timing perturbations can be greatly 
enhanced by many-body dynamics — we have assumed independent Keplerian 
orbits for the subs ystems. Resonant ef f ects on timing perturba tions have 



been considered by [Holman and Murray! [2005l | , [Agol et al.[ |2005 
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Figure 11: Sensitivity of eclipse timing measurements to circumbinary plan- 
ets. The vertical line represents the approximate critical orbit around a 
10-day period binary. The calculations assume the binary consists of two 
stars each massing 1 Mq, 6 hour echpses, Nobs — 25 observations (150 total 
hours of data), V — 10 magnitudes, and l-cr detections. Prom top to bottom, 
lines show sensitivity for D = 0.1 m on the ground, D = 0.5 m on the ground, 
D = 1.0 m in space (i.e. Kepler), D = 10 m on the ground, and D = 2.5 m 
in space (HST, SOFIA). 
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Figure 12: Sensitivity to circumbinary planets, comparing radial velocity 
and eclipse timing techniques. All calculations assume solar mass stars. RV 
assumes 20 ms~^ precision; both the system velocity and apparent period 
variation observables are included for the sensitivity curve, the latter assumes 
a 5 day binary orbital period. Eclipse timing assumes either aim ground- 
based photometric telescope or 2.5 m space-based telescope (such as HST or 
SOFIA), observing a.V = 10 magnitude system, and 25 observations. 

3.3 Observational Precisions 

Figure [12] compares the sensitivity of RV and eclipse timing to circumbinary 
planets. 
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